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^ ! Abstract 

We propose and analyze a way in which effective multicomponent condensates 
can be created inside high-Q multimode cavities. In contrast to the situation 
C^) ' involving several atomic species or levels, the coupling between the various 

components of the dressed condensates is linear. We predict analytically and 
' numerically confirm the onset of instabilities in the quasiparticle excitation 

m ' 

spectrum. 
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I. INTRODUCTION 



The Bose-Einstein condensation of low density atomic samples provides one with 
a new paradigm in many-body theory, atomic physics, quantum optics, and nonlinear dy- 
namics. Below the critical temperature, condensates are described to an excellent degree of 
accuracy by a scalar nonlinear Schrodinger equation, the Gross-Pitaevskii equation describ- 
ing the dynamics of the condensate wave function [[^J. The elementary excitations of the 
condensate evaluated from a Bogoliubov linearization about the condensate solution are in 
good qualitative agreement with experiments |7|-|il| , and so is the Hartree mean- field energy 
of the system [ |l5| , p^] . 

Nonlinear Schrodinger equations are of course ubiquitous in physics, and have been 
studied in great detail in the past, in situations from fluid dynamics [|17]] to phenomenological 



models of field theories JT8| and to nonlinear optics flnj. They play an important role in 
the study of pattern formation in beam propagation | 20f , and their soliton solutions find 
applications in problems such as light propagation in fibers |PH| . From this work, it is 
known that the dynamical and stability properties of multicomponent nonlinear Schrodinger 
equations can be vastly different from those of their scalar versions, and lead to a wealth of 
new effects pl|-F24|1. It would be of considerable interest to generalize these ideas to the case 
of matter waves and to have multicomponent condensates available. 

While it is not generally possible to create two coexisting condensates inside a trap, 
exceptions are possible, as recently demonstrated in Rubidium experiments by the JILA 
group [p5|| . However, this coexistence relies on a fortuitous coincidence of the scattering 



lengths for the two Zeeman sublevels involved [2q|, a coincidence that cannot be generally 
counted on. There are proposals to optically change the s-wave scattering length of ground 



state atoms f27P , but whether this can be used to produce coexisting condensates remains 
to be seen. The goal of the present paper is to propose and analyze a method by which 
effective multicomponent condensates can be generated inside a high-Q multimode optical 
resonator. The cavity photons dress the condensate, very much like atoms can be dressed 
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by electromagnetic fields [p8[ , and the various dressed condensate states are coupled, e.g. 
via an electric dipole interaction. Hence, the condensate inside the cavity should be thought 
of as a coupled multicomponent system, each component subject to a nonlinear equation, 
and in addition coupled to its neighboring components. In contrast to the situation involv- 
ing two (or more) atomic species or levels, the coupling between the various components 
of the dressed condensate is linear, rather than resulting from collisions and hence non- 
linear. Nonetheless, we submit that this method permits to generate and study "coupled 
condensates" in a controlable — and at least in principle simple — way. 

This paper is organized as follows: Section II defines our model and uses a Hartree varia- 
tional principle to derive coupled nonlinear Schrodinger equations describing the evolution of 
a dressed condensate in a two-mode cavity. Section III specializes to the case where only one 
photon is present inside the resonator, and thus only two dressed condensate components 
are of importance. In that case, the problem reduces to the so-called discrete self-trapping 
equations for a dimer familiar in nonlinear physics. These equations are integrable in free 
space, but not in the trap situation that we consider here. We solve them approximately 
in the Thomas-Fermi approximation and study the spectrum of elementary excitations. We 
predict the onset of instabilities in the system, and compare these analytical results with an 
exact numerical solution of the equations. Finally, Section IV is a summary and conclusion. 

II. PHYSICAL MODEL 
A. Basic theory 

Our model system comprises a Bose-Einstein condensate (BEC) which interacts with two 
counter-propagating modes supported by a high-Q ring cavity. We assume a cavity QED 



configuration |29| and neglect all field modes except the two of interest, so that the electric 



field operator can be written as 



E(r, t) = ie£ p ai e ikz + a 2 e~ ikz e~ luJct + h.c, (1) 
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where e is the unit polarization vector of the light, S p = yhu> c /2eoV is the electric field 
per photon for light of frequency u c in a mode volume V, k = u c /c = 2n/X c the light 
wavevector, and a\ 2 , and annihilation and creation operators of the cavity modes 

satisfying Bose commutation relations [a;, a]] = For this treatment we neglect the 
detailed mode structure of the field in the transverse plane perpendicular to the optical-axis 
z, assuming that it is homogeneous on the spatial scale of the BEC, and that it is unaffected 
by its presence. 

The atoms comprising the condensate are confined by an external trapping potential U (r) 
which binds the atoms on a sub-wavelength scale along the longitudinal axis z. In addition 
they interact with the cavity field, which induces transitions between the ground and excited 
electronic states. The single-particle Hamiltonian H Q for the atoms, in an interaction picture 
with the optical frequency removed, then reads [^C 

H = -J— + U{t) + h6a + a- + ^0 Q [(7 + (ai + a 2 ) + (7_(a| + a\)}, (2) 
2m 

where we have located the BEC at z = without loss of generality, p is the center-of-mass 
atomic momentum, m the atomic mass, 5 = u a — uo c is atom-field detuning, uo a being the 
atomic transition frequency, Qq = d£ p /h is the strength of the atom- field coupling, d being 
the atomic dipole-matrix element, and a + , o~_ are pseudo-spin atomic raising and lowering 
operators for transitions between the ground and excited atomic states. We consider the 
case of large atom-field detuning for which the excited atomic state can be adiabatically 
eliminated. This results in the following effective single-particle Hamiltonian (see Appendix 
A for details) involving only the ground atomic state 

p2 ^2 

H e ff = h U(r) -\ (a\ai + ata 2 + a\ai + a[a 2 ) ■ (3) 

2m o v ' 

The four terms involving field mode operators in this effective Hamiltonian describe virtual 
transitions involving the absorption of a photon from mode 1 followed by re-emission into 
mode 1, the same but for mode 2, absorption of a photon from mode 1 followed by re-emission 
into mode 2, and vice versa. The last two of these processes are allowed since the length 5z 



of the BEC is taken to be less than an optical wavelength, yielding an associated momentum 
uncertainty Sp z ~ h/Sz > 2hk. Hence the momentum deficit involved in the transfer of a 
photon from one direction to the other around the cavity is within the Heisenberg uncertainty 
principle. We also note that the effective Hamiltonian conserves the total number of photons 
n = (a\ai + a 2 a 2 ) which is then a good quantum number. 

Proceeding now to include many-body interactions, the second-quantized Hamiltonian 
describing our system is |§ 

H = J tZltffi(l| J ff 6// |2)*t(i)*(2) + iy"d{^}(l,2|y|3,4)*t(i)*t( 2 )*(3)^(4), (4) 

where I denotes a full set of quantum numbers, and and W(£) are the usual atomic 
field annihilation and creation operators, which for bosonic atoms satisfy the commutation 
relations 

= S(£-t'). (5) 

The two-body potential is in the limit of s-wave scattering || 

V(r 1 ,T 2 )=hV 6(r 1 -T 2 ), (6) 

where Vq = Anha/m measures the strength of the two-body interaction, a being the s-wave 
scattering length. Here we consider a repulsive interaction so that a > 0, Vo > 0. 

B. Coupled-condensate equations 

The second-quantized Hamiltonian for our system conserves both the number of atoms 
and the total number of photons, so we consider a state comprising N atoms and n photons. 
The state of the system can be written in the form 

\*N,n(t))= £ dn... dr N f ni 

t)*t( ri )...^t( rjv )| 0)ni)n2 ) ) (7 ) 

n\+n2=n 

where the summation runs over all positive integers ni ( 2 obeying n\ + n 2 = n, 
/m,n 2 ( r i5 • • • j t n, t) is the many-particle Schrodinger wave function for the BEC given there 
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are n\ photons in mode 1 and n 2 photons in mode 2, and |0, ni,n 2 ) is the state with no 
ground state atoms present, n\ photons in mode 1 and n 2 photons in mode 2. To proceed 
we invoke the Hartree approximation, which is appropriate for a Bose condensed system in 
which the atoms are predominantly in the same state. The Hartree approximation is there- 
fore strictly valid at zero temperature, and for a weakly interacting Bose gas as assumed 
here, so that the condensate fraction is close to unity ||. Accordingly the many-particle 
wave function is written as a product of Hartree wave functions 

JV 

/ni,n 2 ( r l; • • • j r N, t) = ]J ni>n2 (rj, t) . (8) 

i=l 

Here 0„ ljn2 (r, t) is the Hartree wave function which represents the state the atoms occupy. 
The equation of motion for ni ri2 (r,t) results from the Hartree variational principle 
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(* N , n (t)\ih%--H\* N , n (t)) 
ot 



0, (9) 



and takes the form of a system of coupled nonlinear Schrodinger equations, or Gross- 
Pitaevskii equations 

h 2 



+ 



V 2 + £/(r) 

2m 



(r, t) + hV N\<p nun2 (r, t)| 2 niin2 (r, t) 

i 

nn 2 



(r,t)), (10) 
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where the photon numbers n 12 in modes 1 and 2 again run over all positive integers obeying 
Hi +n>2 = n. In the limit i7 = there is no coupling between the cavity modes and the BEC 
and Eq. (|10D is the usual scalar Gross-Pitaevskii equation for the condensate. In contrast, 
for non-zero values of Qq the processes involving absorption of a photon from one direction 
and re-emission into the other direction lead to a linear coupling between the state with 
(ni, TI2) photons and those with (rti — 1, n% + 1) and (ri\ + 1, n 2 — 1) photons, the notation 
(rii, n 2 ) meaning n\ photons in state 1 and n 2 photons in state 2. As a result, the system is 
generally a superposition of states with different (rii,n 2 ). 

As they stand, Eqs. fllPp account for the full three-dimensional structure of the BEC. 
In order to make our presentation as straightforward as possible we now make some further 
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simplifying assumptions, but we stress that these are not essential and do not limit the 
generality of the conclusions that we draw. To proceed we write the trapping potential 
explicitly as @-§ 

U(v) = ^1 (ri + A 2 . 2 ) , (11) 

thereby separating the longitudinal potential out from the transverse trapping potential. 
Here r = (r±,z),r± being the transverse position coordinate, u± the transverse angular 
frequency of the trap, and A = uo z /uj±_ is the ratio of the longitudinal to transverse frequencies. 
Here we assume that A < 1 so that the longitudinal trapping is much weaker than the 
transverse trapping, hence giving the BEC density profile a cigar structure fLElfl Specifically, 
we assume that the transverse structure of the BEC is not significantly altered by many-body 
interactions and is determined as the ground state solution of the transverse potential 



hu±v g (r ± ) 

and we express the Hartree wave function as 



h 2 r, muo\ 



-— vl + 

2m 1 



v 9 (*±), (12) 



ni , n2 (r,t) = v g (r ± )e-^ V>' (M)- (13) 



Substituting this expression into Eq. (jlC]), projecting out the transverse mode, and drop- 
ping the prime for simplicity in notation, yields the coupled Gross-Pitaevskii equations 
for the quasi-one-dimensional system. Introducing the dimensionless length £ = z/a 



Zl 



a z = ^h/2rmu z being the characteristic length associated with the longitudinal trapping 
potential, and the dimensionless time r = u z t, these coupled equations can be written in 
the scaled form 



x The opposite case of weak transverse trapping A > 1 corresponds to the BEC having a pancake 
structure [|J and alters only the dimensionality of the resulting equations. In particular it leads to 
a two-dimensional rather than a one-dimensional problem. 
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+ 9 (\lni(n 2 + 1)0 + \/ (ni + l)n 2 </> ni+l,n 2 — 1 (14) 



where the Hamiltonian Hl is given by 

«9 2 I 



Hr 



(15) 



<9£ 2 4 

= fig/5u; 2 is the atom-field interaction energy per photon in units of hou z , which acts as 
the coupling coefficient between different states (ni,n 2 ), and 

NV /V _ NV fdT ± \v g (r ± )\* 
uj z a z u z Jdr±\vg(r ± )\ 2 ' 

is the many-body interaction energy for iV atoms in a volume V in units of %uj z . 

Equations ( |I~4"D are the basis of the remainder of this paper. The transformation to 
dimensionless variables reveals that the key parameters for the system are the linear coupling 
coefficient g, and the nonlinear parameter rj describing self-phase modulation. 

III. DRESSED BOSE-EINSTEIN CONDENSTATES 
A. Dressed BECs 

The dressed BECs are the eigenstates of Eqs. (|14| ) (or more generally Eqs. (|I~0])), and 
are quantum superpositions of states with different photon numbers (ni,n2). Setting 

<P ni ,n 2 (Z,T) = e~^9 nun2 (0, (17) 
for the dressed states, we obtain 

A^m,n 2 (0 = H L 9 nitTl2 + ^|6'n 1 ,n 2 (0| 2 ^ni,n 2 (0 



+ 9 [yni(n 2 + l)0 m _i )Tl2+ i + y (n x + l)n 2 9 

ni+l,n 2 — 1 (18) 

with /i the chemical potential scaled to hu z . Admissable solutions should also be normalized 
according to 
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d£ E I^WOI 2 ^!. (19) 

3 ni+n 2 =ra 

Numerical calculations are generally required to solve these equations for the dressed states, 
but simple limiting cases can be treated analytically. 

B. Single cavity photon, n = 1 

The essential physics of dressed BECs can be exposed using the simplest case of one 
cavity photon, n = 1, so only the states with (1, 0) and (0, 1) are relevant. Then the coupled 
time-dependent Eqs. ([14]) reduce to 

i0oi(£,T) = H L (j) Q1 {i,T) +S , 0io(£,t) +r?|0oi(£,r)| 2 O i(£,r), 

i0io(e, r) = H L <P W (£, t) + #oi(£, r) + v\<t>io& ^)| 2 <M£, r), (20) 

which yields the following pair of coupled equations for the dressed states 

//MO = #lM£) +sM£) + #oi(£)| 2 M£), 

//MO = #lM0 + sMO + #io(OI 2 MO- (2i) 

These systems of equations are similar to those that appear in the theory of multi-component 



condensates |24] . However, instead of a nonlinear coupling due to cross-phase modulation we 
have here linear coupling due to the exchange of photons between cavity modes via virtual 
atomic transitions. In this respect our equations more closely resemble those describing the 
linear evanescent coupling of adjacent nonlinear optical fibers |32 



Equations ( P0|) are known in nonlinear physics as the discrete self-trapping equations for 



a dimer (see e.g. |22fl and references therein). The stationary solutions for such a system 



were classified and their stability was studied in Ref. |22[ in the case Hl = 0. Three types of 
solution for the dimer were uncovered: an in-phase solution denoted (ff) for which, in our 
notation, 8 W = O i, an out-of-phase solution (f|) with W = —9oi, and asymmetric solutions 
(T •) and (• t) with | #io | + |M- 
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The same classification scheme can be employed here for the dressed states including the 
Hamiltonian Hl in Eqs. In this case, however, exact analytic solutions are not available 

for the dressed states, but approximate solutions can be obtained for the in-phase and out-of- 
phase solutions within the Thomas-Fermi approximation in which the nonlinear interaction 
term dominates over the kinetic energy term ||. In the Thomas- Fermi approximation we 
then obtain the following dressed state solution from Eqs. fl2lD 



d to(0 = ^r 2 )J»T9-\e, (22) 

when the argument of the square root is greater than or equal to zero, and is zero otherwise. 
The top sign in Eq. (p2|) corresponds to the in-phase solution, and the lower one to the 
out-of-phase solution. The normalization of the wave function then leads to the following 
expression for the chemical potential of the two solutions 

^ = ±9 + \e m , (23) 
where £ m = [3///2] 1 / 3 is the longitudinal coordinate at which the Thomas- Fermi solution 



vanishes. Using this expression for the chemical potential in Eq. (|22| ) for the dressed state 
solution we readily find that the profile 

|^ | 2 = ^-^ (24) 

is in fact independent of g and whether it is the in-phase or out-of-phase solution. 

C. Elementary excitations 

The elementary excitations of the system can be found by linearizing Eqs. (|20|) around 
the dressed state solutions 

0io(£, r) = e"^[MO + «io(Oe^ T + ^(0^1, 

0oi(£, r) = e"* T [M0 + «oi(Oe^ T + <i(0^ T ], (25) 
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where and v(£) represent small perturbations around the dressed state with energies 
\i ± u. In zeroth order, substitution of these expressions into Eq. ( pH|) results in the system 
of two equations fl2~T|) for the dressed states. In first-order, it leads to the system of four 
equations for the linearized perturbations and t>(£) 



^wio(C) = 




+ 2#io(£)| 2 


-M 




^woi(C) = 


H L 


+ 2v\8io(0\ 2 


-M 


«oi(0 + ^«io(0 + ^oi(OV(0. 


CJfio(C) = 


H L 


+ 2??|MOI 2 




uio(0+^oi(0+^io(OVo(0. 


LOV oi(^) = 


H L 


+ 2v\Oio(0\ 2 




^oi(0 + ^io(0 + ^oi(OV(0- 



The normal modes of this system of coupled equations are identical to the elementary exci- 
tations determined via the Bogoliubov method in which the Hamiltonian for the linearized 
perturbations is brought into diagonal form using a Bogoliubov transformation . 

It was shown in Ref. p2| that for the case corresponding to Hi = the out-of-phase so- 
lution (||) of the self-trapped equations is always stable, while the in-phase solution (||) is 
stable until it bifurcates at a condition corresponding to ji + = 2g, yielding a stable asymmet- 
ric branch and an unstable in-phase branch. Here we study the influence of the Hamiltonian 
Hl on the stability of the in-phase and out-of-phase dressed-states of the system. 

An exact solution for the normal modes of Eqs. (|26|) is not available, to the best of 
our knowledge. To proceed we therefore use the consequence of the Thomas-Fermi ap- 
proximation ^ m > 1, which means that the profile of the dressed state solution is broad 
compared to the characteristic length scale of the trapping potential. This allows us to 
assume 1$^ n2 (£)| 2 — \^nm 2 (^ = 0)| 2 ~ £,m/^V f° r normal modes localized close to the cen- 
ter of the trapping potential. With this replacement Eqs. (B^) can be conveniently solved 
by expanding the perturbations in terms of eigenfunctions of the linear trapping potential 

%(6 = (" + i/M 

MO = E^(0, MO = £%<?„(£); 

V V 

vio(0 = £c?o?„(0, v 01 (0 = £ <&<&,(£). (27) 

V V 
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Substitution of these expressions into Eqs. ( p6|) gives a system of linear equations for the 
coefficients b u , c v which is straightforward to solve. For the out-of-phase dressed-state the 
spectrum of the elementary excitations is 

" = { ± l/(" + ^5a)HWH4*)(" + 5 + 2 ')}- (28) 

Both branches of normal modes given by Eq. are stable, as they are characterized by 
real values of to. In contrast, the spectrum of the normal modes for the in-phase dressed-state 
is 

- = {±^+i+ie-2 S )(,+i-2 9 ), ± ^ + i + ia)(,+i)}. (2 9 ) 

There are again two branches in the excitation spectrum, one of which, 

" = ± ^(" + 5 + 5&- 2 ») (" + I" 29 ) 

can become unstable, that is, us can assume imaginary values. In particular, the region of 
instability is defined in terms of the index v of the linear oscillator mode as 

2<7 - \Ca < v + 2 < 2 9 for 9 > "Cm! 

" + \<2g for g<\e m - (30) 

A detailed analysis of the eigenmodes corresponding to the elementary excitations reveals 
that for the in-phase dressed state the unstable excitations have normal modes that are it 
out of phase, that is, «io(C) = — M oi(0 and fio(0 = — ^oilO? whereas the system is stable 
against symmetric perturbations, Uio(0 = u oi(0 an d v w (£) = t>oi(0- This means that 
as the instability develops, the density profiles of the (1, 0) and (0, 1) components should 
display modulations which are 7r out of phase. 

From this analysis we can also estimate the index v max for the mode of largest growth 
rate, that is, the most negative value of 

fly) = " 2 =(v + \ + \e m -2g)(v + \- 2g) . 
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The requirement f'{y) = then leads to the condition 



max 




(31) 



where the nearest integer value should be taken for v max . The corresponding growth rate is 
readily found to be Im{uS) = r]\8i (0)\ 2 « £ m /^- 



In this section we present sample numerical simulations of the development of the pre- 
dicted instability for the in-phase dressed state BEC. The aim of these simulations is to 
validate the approximate stability analysis of the previous section, and to put its predictions 
in context using a concrete example. 

We have solved the coupled Gross-Pitaevskii equations (p0| ) using a standard beam prop- 
agation technique for the initial conditions 



where #10 (0 = #oi(0 is the in-phase ground state. In the initial condition (j32|) the parameter 
e < 1 is included to provide a slight deviation from the exact in-phase solution. This 
deviation from the exact ground state of the trapping potential can be viewed as a wave 
packet of the normal modes, which triggers any instability present in the system. 

We numerically generated the in-phase ground state solution by evolving the Thomas- 
Fermi ground state (22), which represents a symmetric perturbation of the system and is 
hence stable, until the density profile reached a steady-state. This evolution of the ini- 
tial Thomas-Fermi solution towards the actual ground state occurs since in our numerical 
scheme an absorber is placed at the spatial grid boundaries to avoid unphysical reflections 
of high spatial frequencies: The absorber removes the high spatial frequencies present in 
the Thomas-Fermi solution leaving behind the actual ground state solution for which the 
absorber has a negligible effect. 



D. Numerical simulation of instability 



0oi (f,o) = vT^MO 



(32) 
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For the numerical simulations presented here we have set r\ = 80, and thus £ m = 4.93 ^> 1, 
so the approximations employed here should be valid. Figure 1 shows the density profile 
|0io(£> t)\ 2 for r = (solid line) and r = 2 (dashed line), and g = so there is no coupling. 
In this case the density profile remains close to the initial profile with no sign of any density 
modulations appearing, meaning that the system is stable against density oscillations. Figure 
2 shows the evolution of the central densities |0io(O,r)| 2 and |0oi(O,t)| 2 for g = 7.29. 
Initially, the densities show modulations resulting from the beating of the normal modes 
excited in the initial state, but this is followed by a region of exponential growth for one 
component and decay for the other. This is where the excitation of largest growth rate is 
expected to be dominant, and the predictions of the linear stability analysis can be tested. 
In particular, for the parameters used here we find v max = 8 from Eq. (|3ip . Figure 3 shows 
the density profiles |0io(£,t = 2)| 2 (dotted line) and |0oi(£; r — 2)| 2 (dashed line), along 
with the initial profile (solid line) for comparison. Here we have plotted the densities for 
t = 2 in the region of exponential growth and density modulations signaling an instability 
are clearly seen. In particular, the density oscillations of the two components are tt out of 
phase as predicted, and the oscillations correspond precisely to those expected for the most 
unstable mode with v max = 8. 

These results clearly show that the linear coupling and associated quantum superposition 
of the system wave function have a large effect: In contrast to the stable BEC shown in 
Fig. 1 for g = 0, the introduction of linear coupling via a single cavity photon is sufficient 
to render the iV-atom condensate spatially unstable. 

The density profiles shown in Fig. 3 are appropriate if we determine which cavity mode 
the photon occupies. For example, |0io(£,t)| 2 is the density given that there is one photon 
in mode 1 and none in mode 2. In contrast, if no determination is made of which mode the 
photon occupies then the density profile is Trf(p atom ) = |0oi(£j r)| 2 + |0io(£, t)| 2 , and this is 
shown in Fig. 4 where the density oscillations remain but with sufficiently reduced contrast. 
Here we see the possibility for a delayed choice experiment with a many-body system: imag- 
ine that the system is left to evolve and then released from the trap after which it falls under 
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gravity and its density profile is measured. If we measure which mode the photon occupies 
as the BEC is dropping we expect the large contrast density oscillations, whereas if we don't 
measure which mode the photon occupies we expect the low contrast density oscillations. In 
addition, the decision whether to measure the cavity photon or not can in principle be made 
after the BEC has left the cavity when the BEC and field no longer interact, thus providing 
a delayed choice experiment. Care should be taken, however, that the time at which the 
measurement is performed is within the linear growth range illustrated in Fig. 2. 



IV. CONCLUSION AND OUTLOOK 

In this paper, we have introduced the concept of dressed condensates, which permit to 
create a coupled, multicomponent macroscopic quantum system whose dynamics can be 
vastly different from that of a bare condensate. A number of immediate extensions of the 
ideas presented here can readily be envisioned. For example, one can easily imagine ways to 
create three-component systems, entangled condensates, etc. Such systems will allow one to 
extend many of the ideas related to measurement theory that have been developed in recent 
years in quantum optics to truly macroscopic quantum systems. 
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APPENDIX A: FAR OFF-RESONANCE SINGLE- ATOM 
EFFECTIVE HAMILTONIAN 



In Section I we introduced the Hamiltonian (0) which describes dynamics of a single- 
atom. This Hamiltonian conserves number of excitations in the system, and thus its state 
can be represented as a linear combination of the states with one excitation 

|V(*)> = a o(t)|eOO) + a 10 (t)\glO) + a ol (t)\g01). (Al) 

Equations for the coefficients oti{t) follow from the Schrodinger equation for the state 
and read (here we omit kinetic energy and trapping potential terms) 

iaoo(t) = Sa 00 (t) + Q (a w (t) + a 01 (t)) 
ia 10 (t) = n ocoo(t) 

ia 0l (t) = n a 00 (t) (A2) 

Assuming 5 ^> 1 the excited atomic state can be adiabatically eliminated 

aoo(t) = -^r(aio(t) + a 01 (t)) 
o 

and the resulting equations for the coefficients aoi(t) and «io(t) read 

ia w (t) = — p(ai (t) + a 01 (t)) 
o 

n 2 

ia 01 (t) = —±(a w (t) + aoi(t)). (A3) 
o 

At this point it is straightforward to see that these equations follow from the Schrodinger 
equation for the state 

\i/>eff(t)) = aio(t)|10) + « i(t)|01) 
if the effective Hamiltonian for the system is (S). 
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FIGURES 

FIG. 1. Evolution of the density profile |</>oi(£, t)\ 2 at r = 2 (dashed line) from the initial 
ground state solution (solid line) without linear coupling g = 0. In this figure and all others we set 
r) = 80. 

FIG. 2. Time dependence of the densities |0oi(€? 7 ")| 2 (dashed line) and |</>io(£, t)| 2 (solid line) 
at the center of the trapping potential showing the region of exponential growth for the case with 
linear coupling g = 7.29. 

FIG. 3. Oscillatory pattern developed on top of the ground state profile in the exponential 
growth region r = 2: Ground state solution (solid line), |</>io(£, t)\ 2 (dotted line), and |<foi(£> r )| 2 
(dashed line). 

FIG. 4. Total density |^io(£, r)| 2 + |<Aoi(£,t)| 2 at t = 2 in the region of exponential growth 
for the case that no determination is made of which mode the photon occupies. 
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